Low-energy Antiproton Interaction with Helium 



W. R. Gibbs 

0\ • New Mexico State University 



Las Cmces NM, 88003 



^ : February 9, 2008 

Abstract 

. An ab initio potential for the interaction of the neutral helium atom with antipro- 

I tons and protons is calculated using the Born-Oppenheimer approximation. Using this 

' potential, the annihilation cross section for antiprotons in the energy range 0.01 /xeV 

^ ■ to 1 eV is calculated. 



1 Introduction 



Recent work on the trapping of low-energy antiprotons [J| presented the measurement 
^ I of the annihilation rate on helium atoms which remained in the partial vacuum. It 

C ' was found that the rate increased sharply as the cooling took place and then suddenly 

dropped to the background level. There has been some speculation concerning the 
^ I cause of this sudden decrease in the rate. One possibility could be that there is a 

^ ' barrier of a few millivolts which would only be important when the energy of the 

antiprotons reaches this level. 

While these observations formed the initial motivation for the work presented here 
there are additional reasons why the knowledge of this potential is interesting. 

First, the experiment mentioned above points out the possibility of carrying out 
experiments in which this annihilation cross section could be measured and, in a more 
advanced version, the x-ray transitions might be measured as well. 

Second, there are more general reasons for knowing this potential. For example in 
the particle remnants of the big bang there could well be some antiprotons left in space. 
To the extent that these antiprotons remained free they would have thermalized with 
the ambient background of photons and so would have an average energy corresponding 
to a temperature of around 2.7 K. Their interactions would be with naturally occurring 
hydrogen and helium. 

As a third reason, it is conceivable that there could exist pockets of potential or 
perhaps changes in slope in the antiproton-atom potential which could allow metastable 
states to be formed. These might be associated with configurations of the electrons 
which are strongly deformed from a spherical distribution. Such pockets of potential 
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are perhaps more probable in larger aggregates of matter but this simple case of helium 
provides an entry into the subject. 

There exist some calculations [Q, ^] of a nature related to this reaction but the 
author knows of no calculation of this precise process. 

The calculation of the annihilation cross section as a function of energy requires the 
knowledge of the electrostatic potential between the p and the neutral helium atom, 
the annihilation potential of the nucleus and the solution of the Schrodinger equation 
with these two potentials. 

It is natural to approach the problem of antiproton-helium electrostatic interaction, 
considering at the same time the proton-helium potential. In this way a check on the 
calculation is provided. The calculation of the proton-helium ab initio (variational) 
potential by Kolos and Peek[^ is considered to be the standard reference. 

While one might at first think that the calculations are almost the same (just the 
charge of the external particle is changed), because of the art of choosing the variational 
wave function they can be very different. For instance, Kolos and Peek (following earlier 
work by Kolos and Wolniewicz[^]) chose a trial electron wave function which was 
elliptical in form, thus naturally tending to surround the two positive charges. This 
function is quite appropriate for the two nuclei separated by small or intermediate 
distances. For large distances, where the electrons are nearly spherically distributed 
around the helium nucleus, this wave function can be expected to less efficient. For 
the case of the antiproton this form is perhaps less appropriate. 

Since it is desirable to do both calculations on the same footing, a single trial func- 
tion has been use for both cases. It is based on an expansion in spherical coordinates 
about the helium nucleus. This function is very efficient for small or large distances 
between the two nuclei. It has more difficulty representing the wave function well at 
intermediate distances where it may require a large number of terms. 

In the following section the expressions for the calculation of the variational ratio 
are derived. In section |2| the results of the potentials are given and in section ^ the 
annihilation cross section is computed. 

2 Calculation of the potentials 

The trial wave function chosen for this work is based on a spherical expansion about 
the helium nucleus in terms of spherical harmonics and orthogonal polynomials in the 
(scaled) radial distance of the electrons from the center. It is taken of the same form 
for the proton and antiproton problem so that a direct comparison can be made. This 
partial wave expansion goes over into the polarizability expansion naturally. 

2.1 Variational Integrals 

A completely general form of the wave function for two electrons can be expressed as 
V^(ri,r2) = ^V',,,„,,,,,^,,„,,„,y™nri)l^r*(r2)ini(yi)i:n,(y2)e-5(s'i+2'2) (1) 
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(2) 

where Ln{r) is the Laguerre polynomial of order 2 (usual notation Ln (r)) and yi = 
y2 = T^jo.^ 1 "a" being a variational parameter which sets the scale of the system. 
The quantities '!/'£i,mi,£2,m2,ni,n2 the components of the wave function in this basis 
and embody (reduced by one index, see below) the other variational parameters of the 
calculation. It has been assumed, as usual Q, that the singlet spin state dominates 
the lowest energy configuration leaving the spatial wave function symmetric. 

Since the two-electron wave function can depend only on the relative value of the ^ 
angles, because of the symmetry around the ^He-p axis, the most general form reduces 
to 

^(ri,r2) = 5: VmA/2,ni,n2>^,T(ri)rr(r2)^n,(yi)Ln2(y2)e-^(^i+?^2). (2) 

It is this expression which will be used in the following work. The sums on I and 
n are taken to imax and rimax- The condition that the spatial electron wave function 
is symmetric requires that the V's are symmetric under the interchange (^i,ni) 

The normalization of the wave function is given by 

^(ni + l)(ni + 2)(n2 + l)(n2 + 2)^^,,^,,,,„^,„,. (3) 
The hamiltonian for the problem is 

^2 /i2 2e2 2e2 , , 

H = V? Vi + 1 r ± — ± — , (4 

2m 2m ri r2 |ri— r2| |ri — R| |r2 — R| 

where R is the vector separating the antiproton (proton) from the helium nucleus. 

In order to carry out the variational calculation the computation of the expectation 
value of the trial wave function of each of these terms is needed. 

The radial derivative part of the kinetic energy can be expressed as 

(5) 

where K is an integer. 

Twice the contribution to the kinetic energy for one of the electrons is given by 

(fi _|_ iVn + 2) 

KK{n', n) = 5n,n' ^ 2na„/^„ + 2(n + 2)an',n-i - 2(n + l)6n',n (6) 

so that the contribution to both will be 

K{ni,n'^,n2,n'2) = ^ KK{n[,ni){n2 + l)(n2 + 2)(5„/_ 

+KK{n'2,n2)ini + l)(ni + 2)(5„/ 1 , (7) 



with 







e-yLr,,{y)KMdy (8) 



(n + l)(n + 2)(2n + 3) _ _ , m(n + l)(n + 2 
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bni,n2= ye~^Lni{y)Ln2{y)dy 
Jo 



(n + l)(n + 2) 



n = min of (ni,n2) (10) 



The angular momentum part contributes 



X \£i{ei + l){n2 + l)(n2 + 2)a„j_„/ (5„/^^„2 + ^2(^2 + l)(ni + l)(ni + 2)a„2,n^5ni,ni • 

(11) 

The contribution to the potential energy of the electron-He interaction is given by 



(ni + l)(ni + 2)6„2^^(5„^y^ + (n2 + l)(n2 + 2)6„^y^5„2,n'2 • ^^^^ 



The expectation value of the electron-electron interaction involves calculating 

- / dyidy2y,f (yi)y,7^ (y2)i:„;(yi)i^n^(y2)r--^>^r(yi)^/r(y2)i^ni(yi)i^n;(yi)- 

(13) 



/(2£i + l) (2^2 + 1) ^0,0,0 ^0,0,0 ^m,M,m' ^m,M,m' J, i i /-..x 

7 2^ y (2£i + 1)(24 + l)'"^i.-C'/i ^2,L,^^'"€i,L/, W2,L,^^ J(ni,ni,n2,n2,i.j, (14j 



where 



I{ni,n[,n2,n'2,L) = [ yldyiyldy2e ^y^^y^^i L^:^{yi)L^i^{y2)-^i^Ln^{yi)Ln^{y2) (15) 



The coefficient for the expansion of a product of Laguerre polynomials 

Lni{.y)Ln2{y) = y^^ni,n2,n3-^n3 iv) 

can be found by using the explicit expression of the Laguerre polynomials 

m=0 \ / m=0 

{m + l)(n3 + 2)An,,n2,n3 = (^1 + 2)!(n2 + 2)!(n3 + 2)! 



(16) 



(17) 



ni, 712, "3 



E 



(-1) 



mi+rn2+Tn:i . 



mi + 777-2 + ma + 2)! 



mi, m2, m3: 



(771 - 777l)!(777l + 2)!r77l!(r72 — 7772)1(7772 + 2)!7772!(773 - 7773)!(7773 + 2)!7773! 

(18) 
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With the definition 



Jo Jo 



Ik,n 



, ,poo {k+m)\ 



A; < 
A; > 



we have 



I{ni,n[,n2,n2,L) = 

^^n[,ni,jiAi'2,n2,j2(^ji,ki^j2,k2 [Iki+L-l,k2-L+l + -ffc2+-L-l,fci-L+l] • 

To evaluate the interaction energy with the antiproton (proton) we need 



/•oo a,np-y 



r /•! /"cxD 



where Y is the distance from the He nucleus to the antiproton in units of a alo 
z-axis. 

If we define 



/•oo 

am{Y) = J r'e-^^dt; m > 0, 



then 



ao{Y) = — ; YamiY) = e"^ + ma„_i(F). 



Y 



Note 



m! 



/o 

If m < then the exponential integral is needed 



ara{Y) = 7^(y) 



El (y ) must be calculated to high accuracy. 
With the definition 



PmiY) 



am{Y) m > 
E_miY) m<0 



we can write 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

R (A) 



Figure 1: Proton helium potential. The solid line is from Ref. [§. The solid boxes are from 
the (4,4) calculation, the open circles from the Fade extension and the crosses and x symbols 
are from 0, |^ . 



2.2 Results for the proton potential 

The results for the proton potential are shown in Fig. |l|. Calculations were made 
at various radii for the pair {imax,nmax) = (2,2), (3,3) and (4,4). The (4,4) results 
(solid squares) are shown in the figure. Also shown are the results of ref. H, ^ and the 
parameterization of their potential by Bosanac and Knesaurek Q. The (4,4) results 
give an adequate representation of the previous results except in the minimum where 
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O'.O 0.5 LO L5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

R (A) 

Figure 2: Potential energy of the two-electron system. 

there is a significant cancellation between the electronic potential energy and the direct 
proton- nucleus interaction. 

A [1,1] Pade approximate was used with the three determinations mentioned above 
to estimate the result of the limit of the sequence. This result is shown as the open 
circles in Fig |l[ While this last extrapolation is significant the result gives a satisfactory 
agreement with the previous work. 

A specific comparison was made with ref. § at R=3.704 A. They found -8.32 
meV while the present work gives -8.10 meV for the (4,4) search. 



2.3 Results for the antiproton potential 

The potential energy for the electrons alone is shown in Fig. ^. It is seen to vary 
between the limits Eq at R = (approximately equal to the binding energy of H~) to 
El (the binding energy of isolated helium) at large distances. 

Writing the electron potential as a ratio of two fourth-order polynomials, including 
the direct He-p electrostatic potential, assuming that the potential varies as 1/i?^ at 
large values of R and subtracting the asymptotic value, Ei, the full antiproton-helium 
atom electrostatic potential can be expressed as 

V{R) = ^^^^^ (31) 
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O'.O 0.5 LO L5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

R (A) 



Figure 3: Antiproton He potential. The triangles are the result of the (3,3) calculation and 
the solid circles are from the (4,4) case. The solid line is calculated from Eq. |31]. 

where z = 1/R. 

In the fit Eq = -14.3477 eV and Ei = -78.9847 eV were used. The constant "d" 
(= e^) has the value 28.798 eV . By searching on parameters in this form against 
the (4,4) resuhs, a best fit was found for g = 0.6474 A and /? = 1595 meV-A^. The 
last value can be compared with 1557 meV-A^ from Ref |^ and 1678 meV-A^ from 
Ref @. 

The antiproton potential is shown in Figs. ^ and ^. It is seen that there is no 
apparent change in slope which might aid in the formation of a metastable state. Since 
any such effect would be expected to be small (perhaps in the few meV range) it would 
have to be active at moderately large distances to be visable, say beyond 3 A. We 
can perhaps understand why it is likely that there is no such effect by looking at the 
structure of the electron wave function. 

The variation in the scale parameter "a" with R is shown if Fig. |5[ It is seen that 
the system undergoes a rapid growth for a baryonic separation inside of 1 A. Figure ^ 
shows the first dipole component of the wave function which provides a measure of the 
deformation of the electron cloud from spherical symmetry. It also shows little effect 
outside of 1 A. Thus any possible infiection in the potential curve would occur inside 
of 1 A when the potential is completely dominated by the simple coulomb attraction 
of the He-p system. 
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„30 I ^ ^ ^ ^ ^ ^ ^ ^ ^ 1 

O'.O 0.5 LO L5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

R (A) 

Figure 4: Antiproton He potential. The meaning of the symbols is the same as in Fig. 

3 Annihilation cross section 

The energy of the antiproton at the nuclear surface will be given by the incident energy 
plus the gain due to the acceleration in the coulomb potential. Since (at 1 fm) the 
coulomb energy is 2.88 MeV, in the range of incident energies considered here the 
annihilation takes place at constant energy to a good approximation. The nuclear 
potential was taken to be the product of a purely absorbed single nucleon strength 
with a Woods-Saxon density with radius 2 fm and diffuseness 0.5 fm for 4 nucleons. 

The system was treated from the point of view of a nuclear optical model, i.e. no 
consideration was given to the possibility of the knock-out of electrons or the electro- 
magnetic transition into atomic bound states. In general such corrections might be 
expected to increase slightly the annihilation cross section. 

The problem could be treated in a perturbation approach. In that case one would 
solve the problem of scattering from the purely electrostatic potential generated in 
the previous section and calculate the annihilation rate from the expectation value of 
the potential just introduced. Since the nuclear potential is very short ranged, the 
amplitude will be proportional to the square of the wave function at the origin. It was 
found in calculating this quantity that it had very nearly the same value as the pure 
coulomb case 

Ciirj) = 2^7/(6^^'' -1)-^ ^ -2^r/; r/ = -ac/v (32) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

R (A) 

Figure 5: Variational scale parameter "a" . 



for energies above 10 ^V. From these considerations it is seen that the cross section 
should vary as l/v"^. Since the potential is strongly absorbing the wave function will 
be modified significantly in the region of the nuclear potential so that it is better to 
solve the full equation as now described. 

Since there exists no known analytic form for the solution of the Schrodinger equa- 
tion in these potentials at short range or in the region of the "surface" of the atom, 
the wave function was calculated over the entire region from to 80 A. The step size 
was changed by an order of magnitude 4 times with 20,000 points calculated in each of 
the first four (overlapping) regions and 160,000 points in the region of largest distance 
from the helium nucleus. 

The results for the annihilation cross section are presented in Fig. |^ It is seen that 
there is a change in energy dependence around lOfiV. Above that energy the cross 
section is proportional to 1/u^ as can be expected for a coulomb system. Below that 
energy the cross section varies as 1/v as expected for a low-energy neutral system. 

One can estimate this transition point using the asymptotic form of the potential. 
We can define a radius representing the "surface" of a potential as that value where the 
incident wave number is significantly affected by the potential, say where the potential 
becomes a fraction, f, of the incident kinetic energy. In the present case the radius Rq 
will be defined by 

E= — = fP/K (33) 



2m 
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Figure 6: Dipole coefficient in the wave function. 



The long wave length limit is reached at kRo < 1. Combining these conditions, it is 
seen that the transition should take place around 2/f /i V or about 20 fiV for f=0.1. 
The s-wave cross section above 10"^ eV can be represented by the expression 



0.5367r 



(34) 



Since this is greater than half of the unitarity limit the full cross section cannot be much 
larger than this. While the overall accuracy of the present calculation does not justify 
the three significant figures quoted in Eq. the result is stable to this precision in this 
energy range. Thus if the absolute cross section can be established at any one value of 
the energy it will be known over a wide range. The measurement of the annihilation 
rate on a small, known, amount of helium introduced into a trap could thus provide a 
direct measurement of the temperature over a certain region of energy. 



4 Conclusion 

The authors in Ref. [|| estimate that the pressure in the container was of the order of 
10^1^ Torr. If we assume that the helium and antiprotons are cooled to 4.2 °K then the 
annihilation rate calculated using the cross section from Fig. is ~ 1.2 x 10~^ 1/sec. 
Since a maximum rate of 8 x 10~^ 1/sec was observed, there is a slight discrepancy. 
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Figure 7: Absorption cross section for antiprotons incident on neutral helium atoms. The 
solid curves gives the result due to the sum of all partial waves while the separated dashed 
curve gives the p-wave contribution. The dashed curve almost coincident with the solid 
curve represents the s-wave contribution alone. The short dashed curve gives the product of 
the velocity and the cross section normahzed to the cross section at 1 meV. 



Possible explanations of this discrepancy include the existance of a higher pressure 
than estimated or the presence of heavier material. 

It was seen that no barrier or inflection of the potential occurred. However, the 
possible corrections beyond the Born-Oppenheimer approximation have not yet been 
considered and could influence this conclusion. 

The author gratefully acknowledges discussions with T. Goldman, M. Nieto and M. 
H. Holzscheiter. 
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